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Abstract 

We review several recent developments in the theory and phenomenology 
of deep-inelastic scattering, with particular emphasis on precision tests of 
QCD and progress in the detailed perturbative treatment of structure func- 
tions and parton distributions. We discuss specifically determinations of a s ; 
higher twist contributions to structure functions and renormalons; structure 
functions at small x and resummation of energy logarithms. 
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1 DIS in the collider age 



Deep-Inelastic Scattering (DIS) is the traditional testing ground of perturbative 
QCD. In the last several years, perturbative QCD has become, as an integral part 
of the standard model, a firmly established theory, which allows performing detailed 
and reliable calculations. Current work is focussed on the precise determinations 
of the unknown parameters of the theory, and the development of reliable compu- 
tational techniques. The parameters to be determined include not only the strong 
coupling a s , which is the only free parameter in the QCD Lagrangian, but also all 
quantities which are determined from the nonperturbative low-energy dynamics, 
and thus, even though in principle computable, must in practice be treated as a 
phenomenological input in the perturbative domain. 

In the case of DIS, the low-energy parameters are the parton distributions of 
hadrons. Thanks to the operator-product expansion and the renormalization group, 
or more in general factorization theorems, the coefficient functions which relate 
parton distributions to measurable cross-sections can be computed in perturbation 
theory. Furthermore, the scaling violations of the parton distributions are related 
to the short-distance behavior of local operators in the theory, and can thus also be 
computed reliably^ Deep-inelastic scattering provides therefore the most accurate 
way of determining the parton distributions which are then used as an input in the 
perturbative computation of hard processes relevant for collider physics, as well as, 
through scaling violations, one of the prime ways of determining a s . 

In recent times, the status of the DIS phenomenology has changed substantially, 
largely due to the advent of the HERA lepton-hadron collider, which has enormously 
enlarged the kinematic coverage of the data. In the simplest fully inclusive case, 
the DIS kinematics is fully parametrized by the virtuality Q 2 of the photon which 
mediates the lepton-hadron process, and the total center-of-mass energy 



of the virtual photon-hadron process (which implicitly defines the Bjorken variable 
x). By definition, the DIS region correspond to large Q 2 ; the accessible values of 
x are then limited by the available energy. HERA has thus opened up a large 
previously unaccessible region in the (x, Q 2 ) plane. However, some very interesting 
complementary information has also been obtained recently by fixed-target exper- 
iments, which have provided high-statistics data especially in the polarized case, 
in neutrino scattering experiments, and in charged lepton scattering in the large x 
region, where W 2 is low on the scale of Q 2 . 

The current state of the art can be inferred from Fig. 1, which displays the 
data on the inclusive structure function F 2 (x, Q 2 ) collected by one of the HERA 
collaborations between 1994-1997, as well as earlier fixed-target data. The data 

2 The theory of DIS is a textbook[p] subject. Here we will follow the notation and conventions 
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Figure 1: The proton structure function F%(x, Q 2 ) measured by the HI collaboration 
at HERA as well as by SLACg], BCDMS H and NMCQ fixed-target experiments 
(from Ref.§). 
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are compared to the behavior expected by solving the Altarelli-Parisi evolution 
equations to next-to-leading order (NLO), with a fixed value of the strong coupling 
a s (M z ) = 0.118. Parton distributions are parametrized at a fixed initial scale Q 2 , = 
2 GeV 2 with a dozen free parameters, fitted by comparing to the data. It is apparent 
that despite the enormous magnification of kinematic coverage, there is perfect 
agreement between the data and the QCD prediction, and the agreement does not 
seem to deteriorate even when approaching regions where the NLO computation 
might be expected to be not entirely adequate. In particular, no signs of deviation 
from the NLO QCD prediction are seen in the small x region, which was unexplored 
before HERA. 

It is important to understand what is and what is not remarkable about the 
agreement of data and theory displayed in Figure 1. The fact that it is possible to 
achieve a good description of a large number of data with relatively few parame- 
ters is not especially interesting per se: it may simply indicate that the functional 
dependence of the data on the kinematic variables is sufficiently smooth. 

In fact, there exist parameterizations of F 2 (x,Q 2 ), constructed by experimen- 
tal collaborations for the purpose of data analysis, which interpolate all available 
data [0. These parameterizations do not use perturbative QCD at all: they sim- 
ply fit a functional form to the data, yet they hardly need a larger number of free 
parameters. Furthermore, several recent parameterizations of F 2 at low Q 2 (where 
perturbative evolution is not valid), derived within various non-perturbative low en- 
ergy models f| appear to give a reasonably good description of the F 2 data (though 
generally not quite as accurate as the QCD calculation) even up to values of Q 2 as 
large as 100 GeV 2 , where perturbative QCD should be used instead. Conversely, 
a recent phenomeno logical parameterization ||, inspired by a linearization of the 
perturbative QCD prediction, describes the low x and large Q 2 data with only 4 
parameters, and turns out to work very well even when Q 2 is as low as 0.05 GeV 2 , 
where the perturbative approach is certainly invalid. 

It would be equally unreasonable to take the latter as evidence of the validity of 
perturbation theory at very low Q 2 , or the former as evidence that at large Q 2 per- 
turbative QCD should be replaced by nonperturbative models: what distinguishes a 
perturbative computation such as the one shown in Figure 1 is its predictive power. 
Namely, the data are parametrized as a function of x at a reference scale Qq, and 
increasing the number of parameters only results in a more accurate description of 
the x dependence of the data at that scale. The structure function at any other 
scale is then predicted, by solving the evolution equations. This is to be contrasted 
with the situation in a fit, whether or not inspired by a theoretical model. There, 
no prediction is available: the x and Q 2 dependence are fitted in a certain kinematic 
domain, and then extrapolated outside it. 

The extraordinary success of the perturbative prediction displayed in Fig. 1 
reflects the fact that the meaning of a "test of QCD" has changed in the last decade. 
The correctness and consistency of the theory are no longer at stake: tests of its 
parameter-free predictions (such as the Bjorken sum rule JTO], Sect 2.2, or the small 
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x slope [[□]], Sect. 2. 3) are impressively succesful. The phenomenological emphasis 
is now on precision determination of parameters which must be used as input for 
the reliable calculation of hadron collider processes, and as backgrounds to possible 
new physics. As an example of this, we will discuss in Sect. 2 the determination of 
a s , where the agreement of values extracted from different experiments is now taken 
for granted, and the aim of current work is to obtain an accurate value of a s for use 
in the computation of standard model processes. On the other hand, "testing" the 
theory has taken the meaning of stretching the usual perturbative techniques based 
on hard factorization. This can be done moving towards regions of lower scale, 
where power corrections are important, as we will discuss in Sect. 3 with specific 
reference to DIS at large x. It can also be done by studying regions where more than 
one hard scale is present, as we will discuss in Sect. 4 in the case of small x evolution. 
An altogether different possibility is to stretch the domain of factorization theorem 
to less inclusive phenomena. Specifically, recent data at HERA (and elsewhere) 
have greatly extended our knowledge of lepton-hadron scattering in the diffractive 
region, both within and outside the deep-inelastic domain. A discussion of this 
topic |TIJ is however outside the scope of this review .0 



2 Determinations of a s 

A glance to the current edition of the particle data book|frS| shows that a sizable 
part of our knowledge of a s comes from DIS, and specifically from scaling violations. 
Even though this kind of determination of a s does not have nominally the smallest 
uncertainty, it is, together with the determination from the Z width, the most 
reliable from a theoretical viewpoint. The main recent progress here, in comparison 
for instance to the previous edition [16[ of the particle data book, consists of the fact 
that, even though the experimental uncertainty has been only marginally reduced, 
the central value of a s (M z ) is now in significantly better agreement with the world 
average, and with the value obtained from LEP data. This is partly due to the 
fact that some pieces of data had experimental problems which have now been 
corrected, partly due to the fact that some sources of theoretical uncertainty hadn't 
been appreciated fully, and partly to the availability of new data. Here, we will 
review recent advances in the extraction of a s from inclusive structure function 
data. Of course, a s can also be extracted from jet rates and event shapes. These 
determinations, which are only possible at large center-of-mass energy, in the lepton- 
hadron case have first become possible at HERA; they are however still in their 
infancy and not yet competitive with inclusive determinations p"? 



3 Diffraction in the DIS region can be included in a perturbative framework by introducing 
fracture functions JT^, which allow a separation of the nonperturbative information from the per- 
turbative onep3 and may prove a useful way to arrive at a reliable theory of diffraction. 
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2. 1 Scaling violations of structure functions 



The scaling violations of the structure functions F 2 and F3 (unpolarized) and g\ 
(polarized) are driven by the value of a s according to 

^F 2 (iV,Q 2 ) = ^L s (iV)F— g '+ 7w (iV)F 2 smg ' 



dt x ' 2tt 
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+ cn f A lqg (N)Ag(N, Q 2 )\ + O(aj), (3) 

where F 2 (N, Q 2 ) = Jq dxx N ~ 1 F 2 (x,Q 2 ) is the iV-th Mellin moment of the func- 
tion F 2 (x, Q 2 ), g and Ag are the unpolarized and polarized gluon distributions, the 
anomalous dimensions "f(N) are the Mellin transforms of the Altarelli-Parisi split- 
ting functions P(x), and the factor c = J2i£i e 1 f° r charged lepton and c = 1 for 
neutrino scattering. In order to extract a s we need thus either to be able to sepa- 
rate out the nonsinglet component of F 2 , or to disentangle the singlet, nonsinglet 
and gluon distribution, or to determine the parity violating structure function F 3 
(which only contributes to neutrino DIS). Knowledge of the next-to- leading order 
corrections to eq. (131) allows then to fix the scale at which a s is evaluated. 



The current |15 ] global average value of a s from unpolarized scaling violations 
is a s (M z ) = 0.117 ± 0.002(exp.) ± 0.004(th.). This average includes results from 
global analyses of electron muonfT9| and neutrino|2(J DIS data, i.e. obtained 



by introducing a parameterization of the singlet, nonsinglet and gluon at a reference 
scale, evolving it up to the values of Q 2 at which data are available, and inverting the 
Mellin transform in order to compare the structure function f(x, Q 2 ) to the data. 
In all analyses included in the above average, computations (and specifically the 
solution of the evolution equations) are performed to next-to- leading order (NLO). 
The current value is to be contrasted with the earlier [IB| value a s (M z ) = 



0.112 ± 0.002(exp.) ± 0.004(th.). The central value has gone up because of a 
re-evaluation of the CCFR result from neutrino scattering, which now gives 



a s (M z ) = 0.119 ± 0.002(exp.) ± 0.004(th.) (the previous central value wasf21 
a s (M z ) = 0.111 ± 0.002(exp.) ± 0.003(th.)). The variation is mostly due to a 
re-calibration of the energy of the neutrino beam along with smaller improvements 
and corrections. Because of the smallness of its error, the CCFR determination 
now dominates the value of a s from scaling violations, and we will thus discuss it 
in some detail. 

The peculiarity of the CCFR determination of a s which justifies its small sta- 
tistical error is the fact that the use of a neutrino beam allows a simultaneous 
determination of the structure functions F 2 and F 3 , and thus a reliable separation 
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of the singlet and nonsinglet components. The range of Q 2 covered by the data is 
reasonably high, 5 < Q 2 < 125 GeV 2 , and an estimate of the error associated to 
power corrections is included. 

A shortcoming of the analysis, however, is the relatively crude form of the parton 
parameterization adopted: no attempt is made to fit the detailed x shape of the 
structure function; this may lead to an underestimate of the systematic error. [pT0f| 
A perhaps more interesting potential source of trouble is related to higher order 
corrections, whose impact was not studied (the corresponding error is assumed to 
be the same as determined in the unrelated analysis of ref . |§ ) . 

In fact, the NNLO anomalous dimensions have been recently determined, at 
least for even integer values of 2 < iV < 10. [p2] Since the large N (see Sect. 3) 
and small N (see Sect. 4) behavior of the anomalous dimensions is also known, 
it is possible to reconstruct approximately the full anomalous dimension by inter- 
polation. A NNLO analysis can then be attempted. The CCFR F 3 (i.e. pure 
nonsinglet) data have been recently reanalyzed in this way. |24| The inclu- 
sion of NNLO terms results in a decrease of the (by about 5 units with about 
10 d.o.f.), and in a lower central value of a s by 5a s (M z ) = 0.004. Since this is 
equal to the total theoretical error quoted in ref.[]20|], it suggests that this error 
may actually be underestimated; a careful inclusion of other sources of theoretical 
uncertainty (specifically related to heavy quark thresholds) suggests that a more 
realistic estimate might be A[th., NL0]a s (M 2 ) = 0.006, which at NNLO is reduced 
to A[th., NNL0]a s (M 2 ) = 0.003 (these numbers apply to the pure nonsinglet anal- 
ysis of ref. |2B]]). An alternative NNLO determination from "world" data finds 
a s (M z ) = 0.1163 ± 0.0016(exp.) ± 0.0016(th.). Here, the smaller statistical er- 
ror is due to using a wider dataset, while the smaller theory error comes from a 
less conservative treatment of theoretical errors: a less general parameterization of 
higher twists is used, and higher orders are estimated by different linearization of 
the Altarelli-Parisi solution rather than by scale variation. 

Another important potential source of error is related to the inclusion of contri- 
butions to the evolution equations beyond logarithmic accuracy, i.e. corrections to 
eq. (|3|) suppressed by powers of Q 2 (higher twist corrections). We will discuss this is- 
sue in detail in Sect. 3. In ref.[^(J an estimate of higher twist corrections is included 
on the basis of a renormalon model (see Sect. 3.1), and assigned a 100% error, which 
gives an uncertainty A(HT, NL0)a s (M 2 ) = 0.001. A (much) more conservative es- 
timate is obtained by parameterizing the higher twist corrections by a function 
which is left as an extra free parameter for each x data bin in the fit. In such case 
the theory error is enormously increased to|23| A(th. + HT, NLO)a s (M 2 ) = 0.010, 
which can only be reduced by performing a NNLO analysis, in which case one gets 
A(th. + HT, NNLO)a s (M z ) = 0.004. This error estimate is perhaps a exceedingly 
conservative, since a very large number of extra free parameters is being introduced 
(16 extra parameters). This inevitably inflates the uncertainty unless the statistics 
is very high; notice that this procedure was adopted in the 'classic' BCDMS deter- 



mination [18] of a s . However, there clearly is a correlation between higher twists, 
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higher perturbative orders and value of a s (we shall come back to this in Sect. 3.2), 
which suggests that higher corrections are not negligible to present accuracy. 

The correlation between higher twist and higher orders raises an interesting 
problem: [26] so far, different sources of uncertainty have always been assumed to 
be uncorrelated, and added in quadrature. This assumption however cannot be 
completely correct, either experimentally or theoretically. On the experimental side, 
many sources of systematics are either common to all data in a given experiment 
or highly correlated point by point, and simply determining the error on each data 
point by adding statistical and systematic errors cannot be right. The neglect of 
correlations not only leads to an incorrect estimate of the total x 2 (which indeed is 
often found to be unrealistically mall, i.e. significantly smaller than one per degree 
of freedom), but more seriously it could bias the best-fit value of a s . 

Indeed, a recent reanalysis|2(J of the BCDMS data shows that the central value 
shifts from a s (M z ) = 0.113[TS|] to a s (M z ) = 0.118 if experimental correlations are 
included. Moreover, on the theoretical side it's clear that not all the uncertainties 
discussed above are uncorrelated: for instance (see also Sect. 3.2) if higher twist and 
higher order contributions are not entirely independent, the associate uncertainties 
should not be simply added in quadrature. Rather, if the higher twists are treated 
as free parameters they should simply re-fitted as the scale is varied. Proceeding in 
this way the result improves rapidly with the increase in statistics. For instance, the 
large error found above remains essentially unchanged if one restricts the analysis 



to the CCFR F 3 data alone |4|, but it decreases to A(th. + HT, NL0)a s (M 2 ) = 
0.004 when the CCFR F 2 data are also included Likewise, the error on the 
BCDMS/SLAC NLO extraction decreases from|l§ A(total)a s (M z ) = 0.003 to[|§ 
A(total)a s (M z ) = 0.0017. 

In conclusion, the value of a s from unpolarized scaling violations can be ex- 
tracted from current data with excellent statistical accuracy: the dominant source 
of error is theoretical, related to higher order corrections. Progress is now related 
to a satisfactory treatment of theoretical errors and their correlations, and reliable 
NNLO calculations. 

A determination of a s from scaling violations of the polarized structure func- 
tion gi(x,Q 2 ) has also been presented recently flQ]| . The result is a s (M z 



0.120lQ;QQ5(exp.) t^oo^th.). The statistical accuracy of this result (which includes 
a detailed study of the error due to the way parton distributions are parametrized) 
is justified by the recent remarkable improvement in statistical accuracy and kine- 
matic coverage of the polarized data. []28 Indeed, the error is already dominated by 
the theoretical one: this is however due to the fact that the available data are taken 
at relatively low Q 2 (the most accurate data have an average Q 2 of a few GeV 2 ), 
so that higher order and higher twist corrections are very large. A competitive de- 
termination of a s would require the availability of polarized data at higher energy 
(such as could be acquired at a polarized HERAp^). 

A byproduct of these determinations of a s is the determination of the polarized 
and unpolarized gluon distribution, which is made possible by the availability of 
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structure function data at several scales Q 2 : according to eq. (|3|), knowledge of the 
nonsinglet and singlet structure function at two scales at least determines simul- 
taneously a s and the gluon. Note that the increase in uncertainty on a s due to 
the need to determine the gluon simultaneously is by construction included in the 
statistical error, and thus appears to be subdominant. Scaling violations provide 
the most precise way to determine the unpolarized gluon for x & 0.3 (where it is 
known to about 10% accuracy at scales of a few GeV) , fl30l and the only available 
handle on the polarized gluon (which is only known for values of N ~ 1 and then 
with an uncertainty greater than 50%) [0, |3l| 



2.2 Sum rules 



Rather than using the full set of structure function data in the (x, Q 2 ) plane to infer 
the value of the strong coupling (and parton distributions) from a global fit to the 
scaling violations of all Mellin moments, it may be advantageous to concentrate on 
one single moment, specifically the first moment of F 3 and Q\. The reason is that 
these first moments are proportional to the nucleon expectation value of a conserved 
operator, the vector and axial currents, respectively, with a coefficient determined 
by the operator-product expansion. 

The fact that the operator is conserved means that its expectation value is 
scale- independent (the anomalous dimension vanishes to all orders), and equal to 
a conserved quantum number: the baryon number in the case of F 3 , and the axial 
charge g& (measured in /3-decay) in the case of g\. The relation between the first 
moments of F3 and g% and the respective charges are known as Gross-Llewellyn- 
Smith and Bjorken sum rules, respectively. The Wilson coefficients that relate the 
charges to the first moments of the structure functions have been determined up to 
N 3 LO (i.e. up to O(af) )- ||32|| It follows that knowledge of the first moment of the 
structure function at two or more scales allows one to check the normalization of 
the sum rule, which is a parameter-free prediction of QCD, and its Q 2 dependence, 
which depends on the value of a s . 



Current tests of Bjorken the sum rule |TI| [33| are in excellent agreement with 
the QCD prediction. It is however more interesting to assume the correctness of 
QCD and thus of the sum rules, and use them to extract a s : in principle, a single 
measurement of the moment at one scale provides then a s at that scale to N 3 LO 
accuracy. It is therefore possible to keep under control and greatly reduce the 
error related to higher order corrections, which is dominant in the determinations 
discussed in Sect. 3.1. 

This favorable state of affairs is spoiled by the difficulty of attaining a reliable 
determination of the first moment of the structure function, i.e. its integral over 
the full x range. Since the invariant energy W 2 diverg GS clS X — > 0, it is obviously 
impossible to measure down to x — 0. However, it has been only recently realized 
that the convergence of the integral at small x is rather slow in the presently ac- 
cessible x region. This slow convergence is due to the fact that both polarized and 
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Figure 2: The structure function F 3 (x, Q 2 ); GLS denotes the area under the curve. 
The shaded area is the contribution to GLS from the small x extrapolation (from 
Ref.[ 



unpolarized nonsinglet quark distributions (and thus, correspondingly, g\ and F3) 
appear to grow as 4^ at small x (in agreement with the unpolarized Regge theory 
expectation, but in disagreement with the polarized prediction of a power falloff of 
gi at small x). 

The importance of the small x extrapolation is highlighted in Fig. 2: it is clear 
that the choice of a specific functional form for the small x extrapolation greatly 
constrains the result for the first moment, and can thus lead to an underestimate of 
the uncertainty due to the small x extrapolation. The associated uncertainty has 
been studied in the polarized case in Ref.fT0[| by explicitly varying the functional 



form of the extrapolation. As a result, the value of a s extracted from the Bjorken 
sum rule acquires a very large uncertainty Aa(M z ) = ILq^- the value includes 
statistical and extrapolation errors, but is dominated by the latter. This large error 
overwhelms any advantage which may be obtained from the accurate perturbative 
knowledge of the Wilson coefficient, and makes the extraction of a s from it phe- 
nomenologically not viable with present data. In the unpolarized case a much more 
optimistic estimate (Aa(M z ) = t^ooy) was obtained |34]] by assuming a power-like 
behavior and varying the exponent of the power. In the unpolarized case, however, 
the issue is largely academic, since the need to use neutrino data to determine F 3 
introduce several other sources of large systematic error, specifically related to the 
normalization of the cross section and the estimate of the charm contributions. 

In conclusion, the extraction of a s from DIS sum rules is theoretically very 
clean and could be a nice laboratory to study the convergence of the perturbative 
expansion of the leading twist Wilson coefficients and its relation to higher twist 
terms. p3|| It is however not viable from a phenomenological point of view at present, 



mostly due to the insufficient kinematic coverage of the data which does not allow 
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a reliable determination of the first moment of the structure functions. 



2.3 Scaling violations at small x 

Perhaps the most surprising experimental discovery at HERA is that the very 
strong scaling violations predicted by perturbative QCD at small x are actually 
seen in the data, rather than being cut off by some nonperturbative mechanism: 
in the HERA kinematic range F 2 (x, Q 2 ) rises at small x in accordance with the 
QCD prediction, |35| rather than displaying the set-in of a nonperturbatively gen- 



erated Regge-like behavior, which many expected and pre-HERA data seemed to 
suggest. [j36fl It is natural to think that if scaling violations are strong, then it must 
be possible to determine a s accurately. 

In fact, this turns out to be case also because of more subtle theoretical reasons. 
Indeed, the small x rise of structure functions in perturbation theory is due to the 
fact that the gluon anomalous dimension has a pole in the N plane at N = 1. 
Solving the evolution equations it is easy to see that this pole translates into a rise 
of the gluon distribution and thus the structure function ||35|| F 2 ~ a s xg(x, Q 2 ) ~ 
a s (Q 2 ) exp[27v/^C] (with £ = ln±, ( = In a ^ Q2) and 7 = 12/3 ), stronger than any 

power of In ^ but weaker than any power of x (poles at smaller values of Re N give 
contributions suppressed by powers of x at small x). 

Expanding the anomalous dimensions about the pole leads to a systematic ex- 
pansion of the small x structure function, where the k-th subleading order is sup- 
pressed by exp[C/£l fc,/2 .||37|| The relevant point here is that truncating this expansion 



leads to an evolution equation (which is the small x approximation to the Altarelli- 
Parisi equation) which is local in x and Q 2 ; in particular, at leading order the 
gluon evolves according to a simple wave equation with a negative mass, while the 
(singlet) quark q(x, Q 2 ) is determined by the gluon: 



d 2 [xg(x, Q 2 



l 2 xg(x,Q 2 ) (4) 



' M ; C ' Q)] = ^ ( *,Q2). (5) 

It follows that the leading asymptotic behavior of the structure function at small 
a; is a rise in the product of the two variables £ and (. The scale of the rise is set 
by a s {Q 2 ) (through the definition of Q. The structure function depends only on 
the combination a = an d the asymptotic dependence on a is then universal 
(double asymptotic scaling |37]]). It is important to notice that the asymptotic slope 



of this rise is a universal prediction: since asymptotically the running of a s is 
driven by the leading order (3 function and independent of the value of a s (i.e. 
of A), the asymptotic slope of F 2 is a parameter-free prediction of QCD, and its 
experimental determination provides a basic test of the theory In fact, this is a 
more fundamental test of the QCD dynamics than those provided by the sum rules 
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Figure 3: Double scaling plot of the HERA data published up to 1996. The dashed 
line is the universal (parameter-free) prediction, computed with a s (M z ) = 0.114. 
(From Ref.g|). 



discussed in Sect. 2.2: the value of the slope of the small x rise is fixed entirely by 
the Casimir operator of the gauge group and the number of quark flavors, whereas 
the values of the CCFR and Bjorken sum rule are determined, respectively, by the 
quark baryon number and by the nucleon /?-decay constant. 

However, for the sake of phenomenology it is again more convenient to assume 
the correctness of QCD evolution, and thus of this asymptotic prediction. In the 
region where the asymptotic prediction holds, but two-loop running effects are still 
relevant, it is then possible to extract a s from the observed scaling violations. Such 
a determination appears to be particularly favorable, since these scaling violations 
are strong, and independent of the form of the parton distributions: in particular, 
one does not need an independent determination of the gluon distribution since the 
quark and gluon distribution are asymptotically proportional. There is however a 
major caveat, namely, whether one can trust NLO QCD evolution in the first place. 
This is not obvious, essentially because higher order corrections to the Altarelli- 
Parisi anomalous dimensions are expected to induce a yet stronger rise: the related 
problems will be discussed in Sect. 4. Also, one may think that eventually the 
perturbatively generated rise should be cut off by a nonperturbative mechanism. 

A simple way of making sure that NLO evolution applies is to test for double 
asymptotic scaling (after dividing out universal subasymptotic corrections p8|, p9|), 



to which it reduces at small x. Such a comparison is shown in Figure 3. Notice that 
the asymptotic prediction only depends on the value of a s : a description based on 
NLO evolution thus appears adequate. 

Of course an accurate determination requires comparison to the full NLO pre- 
diction (rather than just the asymptotic behavior). Such a determination was done 
in Ref . |40|, |4l| , on the basis of the data then available, with the result a s (M z ) = 
0.120±0.005(exp.) ±0.009(th.). On top of the advantages already discussed (strong 
scaling violations, little dependence on the input parton distributions) it also turns 
out |Il|] that higher twist corrections are negligibly small in the HERA kinematic 
region (a result which recent data have confirmed ||42||). The relatively large ex- 
perimental error is due to the poor quality of the data then available, which also 
justified a crude treatment of the experimental systematics. The large theoretical 
error is mostly due to uncertainty related to the possibility of deviation from NLO 
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Figure 4: Contribution to the x 2 of the fit from various datasets for the global fits of 
Ref.|43| (left) and Ref.|44| (right, adapted from the original ref.). In the latter case 
the curves, from top to bottom on the left denote the contributions from ZEUS, 
BCDMS, NMC, HI, SLAC, CCFR F 2 , CCFR F 3 and E665. 



evolution at small x, and to factorization and renormalization scale uncertainties: 
these could be both reduced if a s were extracted from the current data (Fig. 1) 
which are much more precise and have a wider kinematic coverage. 

The advantage of an extraction of a s at small x can be seen by comparing the 
contributions to the \ 2 of a global fit of parton distributions to DIS data from 
various datasets as a function of a s (Figure 4).[] Of course, a global fit doesn't 
necessarily provide a good determination of a s for each of the separate datasets 
because of the difficulty in combining the data from different experiments and in 
different kinematic region (and in particular their uncertainties). So, for instance, 
the BCDMS data, which give an excellent ||18|| determination of a s , would appear 



to have no minimum in the fit of Ref. |H4J . Yet it is clear that the minimum in the 
small-x data always is more stable, deeper, and self-consistent than any of the other 
datasets. A determination of a s from the data shown in Fig. 1 is highly desirable 
and could be one of the most precise on the market. 

3 Large x and power corrections 

Large values of x correspond to the kinematical region where the invariant mass of 
the final state W 2 eq. ([!]) is small. In this region power corrections proportional to 
either Aj or ^ (where M is the hadron target mass) are potentially large. Target 
mass corrections are of purely kinematical origin, and are thus completely under- 
stood from a theoretical viewpoint, and routinely included in current unpolarized 
data analyses. Note, however, that in the polarized case the computations neces- 
sary in order to implement the relevant phenomenology have been performed only 
recently. Power corrections in A^ can be generated by the resummation of logs 



of 1 — x which are present at all orders in perturbation theory. 



4 The CTEQ fit@ uses the HERA (i.e. Hl+ZEUS) collected up to 1994 and the old plj C CFR 



data, while the MRST fit|lj includes HERA data up to 1995 and the new CCFR data 
BCDMS j|] data are the same in the two fits. 



201; the 
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3.1 Leading ln(l — x) and power corrections 



Logs of 1 — x are generated in the perturbative computation of the DIS cross section 
for kinematic reasons: as x — > 1 the phase space for the emission of an extra parton 
is suppressed, so that each extra emissions carries a factor of a s (Q 2 ) ln(l — x). At 
the leading logarithmic level, these contributions are resummed|f46H by replacing 



a s (Q 2 ) - a.((l - x)Q 2 ) = Q ff } r (6) 

1 + p a s (Q A ) ln(l - x) 

and bringing the coupling inside the x integration when evaluating the Mellin mo- 
ments of the splitting function on the r.h.s of the evolution equations eq. (|3p. This 
means that the coupling is effectively evaluated at the scale W 2 eq. (0); as x — > 1, 
W 2 — > 1 and the growth of the coupling in the infrared generates a factorial di- 
vergence in the perturbative expansion, which is in turn related to a powerlike 
correction to the structure function. 

This can be understood by considering the, say, first Mellin moment of a s after 
the substitution eq. @: 



r \ oo „i 

I[a s ] = / dxa s ((l - x)Q 2 ) = a s (Q 2 ) ^(-l) fc [/? a,(Q 2 )] fc / dx ln(l - 



x) k 



k=0 

oo 



= a s (Q 2 )J2k\MQ 2 )(3 ) k . (7) 

k=0 

The factorially divergent series can be summed by Borel resummation, i.e. by 
defining 

PCX) 

I[a s ] ee / dse-/ a ' wB) fl[s]. (8) 
J o 

It is easy to see that B[s] = 1 _g /3o - Thus, if we want to determine J[a s ] from eq. (||) 
we must integrate over the pole at s = s = ■h. This generates a contribution to 

I[a s ] proportional to e _s °/ Q!s ^ 2 ' ) = Reconstructing the full x dependence from 
the computation of all Mellin moments shows that the power suppressed contribu- 
tion is actually proportional to i.e. rapidly growing at small x. 

The strength of this 4^ contribution is ambiguous however: for instance, its 
sign depends on whether the path of integration along s is taken above or below 
the pole. The ambiguity is to be expected since we have obtained the result from 
a leading twist computation, i.e. not including terms suppressed by powers of Q 2 , 
and can only be fixed by performing the computation to power accuracy. Of course, 
as long as we are interested in the leading twist, leading ln(l — x) result, we may 
simply refrain from expanding out az s ((l — x)Q 2 ). 

This suggests two directions of progress in the description of the large x region. 
On the one hand, we may push the resummation of ln(l — x) contributions beyond 
the leading log level. Indeed, the next-to-leading resummation is also known; B?] 
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and a systematic method to resum subleading contributions to all orders has been 



suggested recently. [[48[ This determines the large N behavior of the anomalous di- 
mensions and coefficient functions. On the other hand, we may use the fact that 
the all-order summation of leading twist contributions generates power-suppressed 
terms as a means to acquire information on the higher twist terms in the operator- 
product expansion. Considerable progress has been made recently in the latter 
direction by means of the so-called renormalon method.^ 

3.2 Renormalons and phenomenology 

Because the structure function is a physical observable, the ambiguity proportional 
to ^ generated by the Borel resummation of classes of leading twist contributions 
must be cancelled by an equal and opposite ambiguity in the next-to-leading twist 
contribution. Indeed, such an ambiguity arises when regulating the power ultra- 
violet divergence of the next-to-leading twist, i.e., the pole in the Borel transform 
must cancel between leading and next-to- leading twist. We can then take the co- 
efficient of the pole in the leading twist as an estimate of the full next-to-leading 
twist computation. The rationale for this is the same as when estimating a NNLO 
corrections by varying the factorization scale in a NLO computation. Of course, we 
would expect this estimate to be only qualitatively correct. 

A specific class of factorially divergent contribution, related to the small momen- 
tum region of integration on loop momenta (infrared renormalons) can be isolated 
and computed for a variety of observables, and may be used to this purpose. 



In the case of DIS, the renormalon contributions to the coefficient functions which 
relate structure functions to parton distributions have been computed both in the 
nonsinglet |50[] and singlet fl5ij sector (the contributions to the nonsinglet splitting 



functions are also known [5^], but have not been used for phenomenology). One can 
then convolute this ^2 contribution to the coefficient function with parton distri- 
butions taken from a NLO, leading twist fit of structure functions data, and view 
the result as an estimate of the O (^^) P ar t of the structure function itself, up to 
an unknown coefficient. The result can be compared to available phenomenological 
fits,fni where the structure function is parameterized as 



F 2 (x,Q 2 )=F l 2 \x,Q 2 ) + ^ (9) 

and the function h(x) is treated as a free parameter for each x data bin. 

After adjusting the coefficient (to about three times of what would be obtained 
by Borel transform integrating below the pole, as discussed in Sect. 3.1) the agree- 
ment with the data appears to be excellent: in particular, the shape (as a function 
of x) of the observed twist four (i.e. O (^2)) contribution to F 2 appears to agree 

5 It is impossible to cover here the very rich literature on this subject. The reader is referred 
to the excellent recent review Ref . fE9f . 
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Figure 5: Higher twist contribution eq. @ to F2 from the fit of ref.[^]. The curves 
are labelled from top to bottom on the right. 



very well with the renormalon estimate. This poses a problem, since quantitative 
agreement is found, whereas only a qualitatively correct description was expected. 
The result is especially surprising in view of the fact that higher twist contributions 
correspond to operators which measure multiparton correlations in the target and 
should thus be target-dependent, whereas the O (^2) contribution computed from 
renormalons is universal, since it comes from a coefficient function. 

A hint on the possible origin of this state of affairs comes from a recent phe- 
nomenological analysis^] of the BCDMS0, SLAC|] and CCFR|U| F 2 data (see 



Fig. 5). The fitted higher twist contribution agrees well with previous fits,[[L8| even 
though here a different 0] set of parton distributions and a wider dataset are used. 
However, if the value of a s used in the computation of the leading twist part is 
changed from the low value used in Ref.|T8|| to a larger value, a substantial reduc- 
tion of the higher twist contribution is observed. If F2 is computed by folding the 
NLO partons with the NNLO coefficient function ||55|| the higher twist is further re- 
duced. A decrease of the fitted higher twist contribution as more perturbative orders 
are included has also been found in an analysis of the CCFR F 3 data; while the 
strong anticorrelation between the value of a s and the size of the higher twist con- 
tributions has been observed in a recent reanalysis of the BCDMS/SLAC data[K 
and in a phenomenological extraction of higher twists from "world" data. ]56) 



This means that the fit doesn't really distinguish between log and power be- 
havior: the higher twist correction extracted from the fits is in large part a pa- 
rameterization of missing leading twist contributions, because of the truncation to 
NLO, or because the value of a s is too small. This explains naturally the success 
of the renormalon calculation, which is equal to the resummation of higher order 
contributions, and thus provides a good approximation to such missing higher order 
terms. Notice that the leading ln(l — x) terms discussed in Sect. 3.1 are included 
in the renormalon estimate. But this also means that the "dynamical" (i.e. target- 
dependent) part of the higher twist, which is unrelated to higher perturbative orders 
and would spoil this agreement, must be subdominant (i.e. the higher twist is "ul- 
traviolet dominated" ||57||). In principle, it is alternatively possible that dynamical 
higher twists have a similar functional form as the renormalon higher twist, and 
are thus not easily disentangled from them. This is presumably true in the x — > 1 
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Figure 6: Comparison of the F2 data with the NLO calculation based on the par- 
tons of Ref . ||54|| , without power corrections (upper), with target mass corrections 
(middle), with target mass and renormalon corrections (from ref.|53|]). "Empirical 
HT" refers to a model functional form for h(x) eq. ([3J). 



limit, in that a behavior is expected to be generic on kinematic grounds, and 
it indeed happens if dynamical higher twists are e.g. estimated in a simple model 
based on free-field theory. |58| It seems rather unlikely in general though, and it 
appears more natural to simply conclude that dynamical higher twists are smaller, 
and hidden by the renormalon contribution. 

All this has also interesting implications for the leading twist phenomenology: 
it implies that the value of a s extracted from log scaling violations can be under- 
estimated if higher twists are large in the fit and, more importantly, the error on 
a s can be underestimated if higher twists are not allowed to vary enough. These 
problems should be mitigated in a determination of a s at small x, where W 2 is 
large and indeed higher twists appear to be small. |]40|, |42|| 

One is thus led to conclude that there is a hierarchy in power corrections to 
structure functions at large x (Fig.6): target mass corrections are characterized by 
a relatively large scale of order of the target mass M ~ 1 GeV; higher twists due to 
the resummation of higher orders in the MS scheme have a characteristic scale of 
several hundreds of MeV; and "dynamical" higher twist corrections appear at a yet 
smaller scale. Even though there is no understanding of this subdominance, it is 
not necessarily problematic: one would expect the natural scale of dynamical higher 
twists to be of order A, and current data are compatible with this. If however the 
scale turned be much smaller, this would pose an interesting dynamical problem: 
naive perturbation theory would appear to work at low Q 2 better than it ought to. 
Be that as it may, the phenomeno logical extraction |59|] of the dynamical higher twist 
corrections, which would give access to parton correlations in the nucleon, appears 
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to be a very difficult challenge which will require new experimental information. 



4 Structure functions at small x 



As we discussed in Sect. 2.3, the rise of F 2 (x, Q 2 ) at small x observed at HERA is 
in spectacular agreement with the prediction obtained by approximating the NLO 
anomalous dimensions with their highest rightmost singularities in iV-space (leading 
singularities, henceforth). This immediately raises a problem: higher order (in a s ) 
contributions to thr anomalous dimensions contain higher order poles, associated 
to higher order powers of In -. One has thus to face the problem of summing these 

contributions. Because at small x according to eq. (|I|) W 2 ~ — , this is related to 
the determination of the leading high energy behavior of parton-parton scattering 
cross sections. 



4-1 Altarelli-Parisi equations and summation of In - 

While solving the Altarelli-Parisi sums leading logs of Q 2 , the inclusion of suitable 
higher order terms in the anomalous dimension allows one to also sum other logs: 
for instance we discussed in sect. 3.1 how to sum leading logs of 1 — x. The lead- 
ing singularities of the anomalous dimensions in the gluon sector have the generic 



form [60] 



oo 



fi = Y j a k sl {k)ij (10) 



k=0 



n=l 



7 W« = Sa f)«(C^) (ii) 



N- 1 



.ij 



where C = 121og2 and the coefficients satisfy lim^oo = 1. At leading order 
in a s , the only nonvanishing coefficients are in the gluon sector, where they satisfy 
a^ 9q = ^a^ 99 (with n > 1), and can be extracted|Kj from knowledge of the 
leading log high energy behavior of QCD perturbative cross sections, which is in 
turn found by solution of the BFKL equation. |n] If such contributions are included 
in the anomalous dimensions, the solution to the Altarelli-Parisi equation contains 
all leading logs of -, because a k-th order pole in N space corresponds to a In - k 
contribution. 

Naively, one may think that since the series in eq. ( |Tl~]) has unit radius of 
convergence, one must resum all orders in In - whenever the expansion parame- 
ter E(x, Q 2 ) = Ca(Q 2 )\n^ ~ 1. This is problematic because C is very large so 
S(x,Q 2 ) ^> 1 in the HERA region, and in fact S(x,Q 2 ) > 1 also in most of the 
SLAC and all of the NMC region (see Fig. 7): the small x region is outside the ra- 
dius of convergence of the series. This, in the pre-HERA age, led to the conclusion 



that the perturbative expansion eq. ( |TT| ) breaks down at small x. ||62|| 
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This conclusion is however unwarranted, because the quantity of direct physical 
interest is the Altarelli-Parisi splitting function, related by inverse Mellin transfor- 
mation to the anomalous dimension eq. (p3[-|TT|). Now, the inverse Mellin transform 
is essentially the same as the Borel transform discussed in Sect. 3.1, so if the anoma- 
lous dimension as a series in a s eq. (|TT| ) has a finite radius of convergence, the split- 
ting function has an infinite radius of convergence for all finite x. It follows that the 
leading In - summation converges for all x, and can thus be used down to arbitrarily 
small x: in each (x, Q 2 ) region, it is sufficient to sum the first k terms in eq. (|TTD, 
where Q 2 )] k /k\ ~ 1, i.e. (from Stirling formula) k « 2.75 for H *S 1. Notice 
the difference in comparison to the large x case, discussed in Sect. 3.1, where the 
leading large N anomalous dimension diverges factorially, so the leading ln(l — x) 
series has a finite radius of convergence (cfr. eq. ([])), one must resum an infinite 
number of terms, and the resummation breaks down when x — > 1. 

In fact, it can be shown [p^] that the Altarelli-Parisi equation with leading- 



singularity anomalous dimensions 

7 (o),<i eq- 

fllCf) is completely equivalent, up to 
higher twists, to the BFKL equation (in the inclusive case), which describes the 
leading log high energy asymptotics. It follows that one may proceed in two equiv- 
alent ways: either one determines the LO coefficients aPW eq. (O) from the BFKL 



kernel, and then one inverse-Mellin transforms the series eq. (|TT| ) term by term, 
thereby obtaining a convergent series which may be truncated at finite order. [Q 



Or, alternatively, one constructs "nonperturbatively" , i.e. numerically the leading- 
singularity anomalous dimension as a function of a/N directly from the BFKL 
kernel. 



Either way, one is lead to a perturbative expansion of the small-x splitting 
function 

pV(x, Q 2 ) = P (0)lj + a s P {1)ij + ... 

1 OO / 1 1 \ n 

pW^^QV^ft ■ (12) 

In the quark sector, the coefficients start at 0(a s ) and are thus subleading 
compared to the gluon ones; they are however also known. f66[ It is thus possible to 
determine the splitting functions P* J to leading nontrivial order, add them on top 
of usual NLO splitting function (with suitable matching conditions), and compare 
the result to the HERA data|>§ §§. 

Unfortunately, the results of this comparison are discouraging: as long as the 
HERA data were still relatively scarce, no evidence of the log - summation could be 



seen,|>|] but as soon as they became accurate enough, the presence of such effects 
could be excluded to very high accuracy, essentially because they induce much 
stronger scaling violations than observed in the data.[]67| In fact, the only way 



current data can be reconciled with evolution equations that include a summation 
of small x effects is by artificially fine-tuning the factorization scheme in such a 
way as to remove these effects from the measured region. 
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Figure 7: (Left) Contours of constant H = Ca s (Q 2 ) ln(-) = 1, 2 ... 20 (from bottom 
to top) in the (x, Q 2 ) plane, as a function of x and Q 2 . (Right) The ratio of NLO 
to LO splitting functions eq. (11) in various factorization schemes: MS(solid), Qo@ 
(dotted), "physical" |D (dashed and dot-dashed) (from Ref.|7 



A possible way out of this impasse consists of exploiting the wider factorization 
scheme freedom which is related to the summation of small- x contributions: for 
instance, one can always reabsorb any undesirable scaling violations induced by the 
leading-order contributions in eq. (P~D[) , which only appear in the gluon sector, in 
a redefinition of the gluon normalization, at the only cost of changing the next-to- 
leading coefficients af£> 99 . It could then be that the effects which cause disagreement 
with the data are a spurious consequence of the choice of factorization scheme, and 
would disappear in a "physical" scheme where each parton distribution is identified 
with a physical observable. [|68| In order to check whether this is the case, however, 
the full set of NLO coefficients d% eq. (|ll]) is needed: one may then verify whether 
in this scheme the undesirable scaling violations actually disappear. 

Thanks to the recent calculation of the next-to-leading corrections to the high- 
energy QCD asymptotics, ||69|| which allows a determination of the next-to-leading 
coefficients a^ 99 , it is now possible to resolve this issue. The outcome of the calcu- 
lation, however, brought a new surprise: the NLO coefficients are negative and grow 
very large [fT0|, [7I| in comparison to the leading order ones as the order k in eq. ( p~2| ) 



grows. In fact, it can be shown both numerically ||70|| and analytically |72] that the 



ratio of the NL contribution to the splitting function to the leading order grows lin- 
early as a function of S, and in fact overwhelms the leading order for any reasonable 
value of x and Q 2 (fig. 7). The problem is in fact worse in "physical" schemes where 
the ratio grows as S 3 / 2 ; the same happens in the so-called Q scheme |73[ (designed 
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to reduce perturbative corrections in the quark sector) , as well as in the momentum 
conserving scheme |74j], though the slope of the growth is somewhat smaller in the 
latter case. 

The growth of the NLO "correction" with respect to the leading order means 
that the NLO P^ eq. (|12"D, though formally subleading, is actually enhanced at 
large H with respect to P^: even though it is suppressed by an extra power of a s , 
at sufficiently small x it always takes over due to the growth of S with In - at fixed 
a s . Hence, it is the perturbative expansion of the leading In - contributions to the 
Altarelli-Parisi equations defined by eq. (12) which breaks down. 

On the one hand, this is good news, in that it explains why the inclusion of 
the leading In - summation P^ u spoils the agreement with the data: these terms 
are actually not leading in the HERA region, and thus there is no reason why 
including them while not including more important contributions should improve 
the phenomenology. On the other hand, however, the fact that standard NLO works 
at all — not only at HERA, but also in the SLAC or NMC region (see Fig. 7) — is 
now puzzling: there must be very large cancellations in the anomalous dimension 
eq. (|TT|), or, otherwise stated, the expansion eq. fll2]) is pathological. 

It is thus necessary to reorganize the expansion. Now, it can be shown |72| 
that the rise of the NL correction is entirely generic — it could have been predicted 
without waiting the outcome of the Fadin-Lipatov calculation: the ratio of the N fc LO 
to the LO splitting functions eq. (|I2]) behaves as |^Jy ~ [a s ln(^)] 2fe_1 . This behavior 
is related to the fact that generically the anomalous dimensions 7^ eq. ([□]) has a 
(2k-l) — th order pole at = 1. The slope of the fc-th order rise is determined 
by the residue of the corresponding pole, which can in turn be simply expressed in 
terms of the k-th order corrections to high-energy asymptotics. 

It is possible to eliminate this singularity in the anomalous dimension, and thus 
the instability of the perturbative expansion, by suitably choosing the factorization 
scale order by order in a s . This is however a fine-tuning: even small variations of the 
factorization scale around the optimal scale bring back immediately the instability 
(see Fig. 8) A better option consists of reabsorbing into the leading order p(°) all 
contributions which lead to the small x growth of the formally subleading PW. It 
turns out that this can be done in a completely unambiguous way by subtracting 
terms related to the residues of the poles in the anomalous dimensions. 

The overall ambiguity introduced by the procedure is contained in a single pa- 
rameter A, which could only be determined by resumming the splitting function to 
all orders. This parameter governs the all-order asymptotic behavior of the splitting 
function P{x,Q 2 ) eq. P{x,Q 2 ) ~ (ln±) 3 / 2 x A . The ensuing splitting functions 

can be viewed as a resummation of the original series: the leading order resummed 
P^ depends on the all-order resummed parameter A, but the ratios of all higher 
order pW to P*- -* is fixed, and tend asymptotically to a constant as x — > (Fig. 8). 

In an Altarelli-Parisi framework this does not entail loss of predictivity, since 
the Q 2 -independent part of the x dependence of the structure function cannot be 
determined anyway. However, it does show that even assuming that the asymptotic 
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Figure 8: The ratio of NLO to LO splitting functions eq. (11) in the DIS scheme 
after resummation (solid), with fine-tuned scale (dashed) and scales a factor of two 
either side of the fine-tuned one (dotted), (from Ref.|fT2[) 



small x behavior of structure functions is determined by perturbative QCD asymp- 
totics, it cannot be determined in any finite order computation and would require 
an all-order resummation of the small- x expansion. 



4-2 Energy evolution 

The resummation of the the leading In - contributions to structure functions dis- 

° X 

cussed in the previous section is not the end of the story, however. So far, we have 
discussed leading In - contributions which are generated when solving evolution 
equations in Q 2 , i.e. the Altarelli-Parisi equations. This, however, is an indirect 
manifestation of a more fundamental phenomenon, namely, the fact that the high- 
energy behavior of the gluon-gluon scattering cross section in perturbative QCD is 
dominated by energy logs. 

In DIS energy logs appear as lnp^ 2 ~ In-, and can be summed by means of 



an equation for evolution in -. This is accomplished by the BFKL equation, [61] 
which in the fully inclusive case can be written in the same form as the usual 
QCD evolution equation, but with the roles of the two kinematic variables x and 
Q 2 interchanged. To this purpose, it is necessary to take a Mellin transform with 
respect to Q 2 and define kernels x(M) as functions of the associate variable M, 
which can then be used to evolve in £ = In - : 

' X 

^G(x,M) = a sX (M)G(x,M) 
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Figure 9: The kernel x(M) eq. fljjD plotted as a function of M. The curves from 
top to bottom in the middle are the LO result with a s — 0.1 and the NLO result 
with a s = 0.1, 0.05, 0.03, 0.01, 0.005. (adapted from Ref.0) 



G(x,M)=f^(^Pj ^G^Q 2 ) (13) 
X (M) = X {0) (M) + a sX (1) + ..., 

where G(x, Q 2 ) is the leading eigenvector of the anomalous dimension matrix eq. ( |T(]| ) 
(i.e. the one associated to the only nonvanishing leading-order eigenvalue), which 
at leading order is just G(x, Q 2 ) = xg(x, Q 2 ) with g(x, Q 2 ) the gluon distribution.^ 
Since G(x, Q 2 ) behaves as a constant at large Q 2 (Bjorken scaling) and falls lin- 
early as Q 2 — > (for kinematic reasons [l]]), the (leading twist) physical region in 
M space is < M < 1. At leading order, x (0) ( M ) = 2 7 - ip(M) - - M) 
(where ip(M) is the digamma function) — the celebrated BFKL kernel — and the 



recent FL calculation [69] in the inclusive case amounts to a determination of the 
next-to-leading correction x {M). 

The qualitative behavior of the kernel is changed dramatically by the NL cor- 
rection (see Fig. 9). To understand the effect of this change, solve the evolution 
equation eq. ( |i"3"D and construct G(x, Q 2 ) by inverting the Mellin transform: 

The asymptotic behavior of the integral at small x, i.e. large £ is determined by 
saddle point. Since the integration path runs along the imaginary axis, the saddle 
is a minimum of x(M) along the real axis, which at leading order is located at 
M = |, where ^ (1/2) = 41n2. But at NLO there is no real minimum, and the 
real saddle is thus replaced by two complex saddles (which must occur in complex 
conjugate pairs by definition of the kernel as the Mellin transform of a real splitting 
function) . |7D|, [75] Hence the NLO, rather than being a small correction to the LO, 
changes the asymptotic behavior, and in fact leads to an oscillatory, unphysical 

6 Eq. ([[3]) applies at fixed coupling a s . Beyond leading order the running of the coupling a s 
with Q 2 must also be included; this introduces some technical complications in the form of eq. ([l3]) 
and its solution which arc however inessential for our discussion. 
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behavior at large £. In particular, the value of x(l/2) has no special meaning if \ 
is computed at NLO. Of course, at small enough a s i.e. high enough scales, the 
leading order behavior is recovered: however, this only happens (Fig. 9) at values 
of a s ~ 0.01 which correspond to a scale above the Planck mass. |70|] 

The original calculation on which these results are based, |69| which is the re- 
sult of an effort of years, has been in large part checked independently: most of 
the individual real and virtual amplitudes have been recomputed with different 
techniques; [[FG] and the elaborate computation that takes from the individual am- 
plitudes to the kernel has also been checked. [77]] It thus appears unlikely that the 
result |)9| is erroneous. The saddle point estimate of the asymptotic behavior has 



also been checked by more accurate calculations, which have confirmed the unphys- 
ical oscillatory behavior of the cross-section. |78|, [79]] 



The problem posed by this unphysical behavior is separate from that discussed 
in Sect. 4.1. Indeed, there we saw that the naive small x expansion of the anomalous 
dimension must be reorganized in such a way that the terms which contribute to the 
asymptotic small x behavior to all orders in the expansion eq. flTD| ) are resummed 
into a parameter A, and included into the leading order of a resummed expansion. 
Now, we are studying how the all-order asymptotic behavior can be derived from 
X, and we see that if x is computed to NL order, the associate asymptotic behavior 
is unphysical. It therefore looks like a further resummation of x is required. 

Several partial resummations of formally subleading corrections to x have thus 
been suggested. First, one can attempt to resum logs of Q 2 related to the running 
of the coupling[|^J, though this actually seems to make things worse. |7TJ Also, one 
may try to optimize the scale choice by the BLM method|81|: this goes in the right 



direction, but does not completely remove the instability. A further option consists 
of noting that the definition of the leading energy logs which are being resummed 
requires a choice of scale: for instance, in DIS at small x, s = W 2 ~ — , and the 
resummation of In - may be more generally viewed as a resummation of In — with 

X * ■ * Sq 

s o — Q 2 - Different choices of Sq then lead to different subleading corrections, and 
a resummation of these "scale-dependent" corrections may change the properties 
of the evolution equation eq. (|i~3|) . f82fl A specific resummation of such corrections 
has been advocated on the grounds that it removes some undesirable singularities of 
x(M) at M = 0, thus leading to an improved asymptotic behavior^, |83|]. This does 
not, however, lead to an unambiguous prescription. Use of an evolution equation 



based on angular ordering seems to further increase the theoretical ambiguities. |84 
A perhaps more fundamental possibility |85|, [86| is to re-examine the assump- 
tions on the high-energy asymptotics that go into the NL determination of x- m 
particular, by introducing physically motivated kinematic cuts in the course of the 
determination of x(M) one can reshuffle contributions between the NL x an d 
the LO kernel x while only changing the (unknown) NNL terms. Again, the 
shortcoming of the approach is that it introduces an ambiguity, in the form of a 
dependence on a "cutoff" parameter. However, the ensuing NL result is then well- 
behaved in a reasonable range of values of this cutoff. The ambiguities can be 
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somewhat reduced ]87|] by again resumming "scale dependent" corrections along the 



lines of Ref. 82 



In conclusion, the ambiguities encountered when trying to construct an energy 
evolution equation appear to be related to the choice of the appropriate large scale. 
Indeed, on the one hand the large log which is treated as leading and summed by 
the evolution equations is £, associate to the large scale W 2 , but on the other hand 
the hard scale which guarantees the validity of perturbative factorization and the 
smallness of a s is Q 2 . The need to reconcile these two large scales appears to be at 
the origin of the difficulties of perturbation theory at small x. A radical suggestion, 
based on the idea that W 2 should also be used as a factorization scale, was discussed 
in ref.f63|; a proof of the corresponding factorization theorems would however be 
required in order for this approach to be viable beyond leading order. We must 
conclude that no reliable theory of small x evolution is yet available. The success 
of simple NLO perturbation theory and the small x approximation to it still lacks 
a completely fully satisfactory explanation. 



5 The embarrassing success of QCD 

The application of perturbative QCD to inclusive deep-inelastic scattering is em- 
barrassingly successful. In many ways, QCD is enjoying the same sort of success 
as the electroweak sector of the standard model. Whether the lack of evidence for 
nonstandard effects should be considered a triumph or a disappointment is of course 
matter of opinion. However, in many instance things are now too easy: we do not 
see deviations from the simple leading log behavior even in kinematic regions where 
we might have expected them, in particular at very large and very small x. This 
reveals that our understanding of perturbative QCD and its limitations is still not 
satisfactory from a theoretical viewpoint. 
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